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We review the consequences of the sensitivity of the relation between the moments 
of a model for the nucleon quark distribution amplitude and the coefficients of its 
polynomial expansion. Criteria for a simpler approach to constructing a model for 
the quark distribution amplitude are formulated. We describe how such a simpler 
(or haplousterotic) model for the quark distribution amplitude of the nucleon is 
obtained from the QCD sum-rule moments of COZ. 



1 Introduction 

Although the knowledge of the nucleon wave function in terms of its fundamen- 
tal quark and gluon degrees of freedom is of outstanding theoretical interest 
because of its process-independence^ there are many drastically different model 
amplitudes for the nucleon available as shown in Fig. ^ All the functions 
shown are intended to describe the distribution of the longitudinal momentum 
fractions of the quarks inside the nucleon. 

In calculations using perturbative Quantum Chromodynamics (pQCD), 
not the full wave function tpt^ itself is used, but the quark distribution ampli- 
tude (QDA) which is defined from the nucleon wave function by 

(f>N{x,fi'^) -.^ I [d^k±i] ^pN{xi,k±i) . (1) 

From Eq. ^is clear that is also process- independent and can thus be used in 
calculations of various observables stirh as formfactors, decay widths or amplk 
tudes in virtual Compton scattering3 (VCS). In the hard-scattering aproachd 
of pQCD the contributions to such processes are assumed to factorize into a 
process-specific hard-scattering amplitude Th between the distribution ampli- 
tudes of in- and outgoing states. For example, a typical contribution to Th in 
VCS looks like the diagram in Fig. ||. 

Whereas the hard-scattering amplitude can be calculated perturbatively 
by exchanging hard gluons, the QDA is a purely non-perturbative object 
and must be obtained from elsewhere. 

Usually, the ansatz for 0jv was an expansion into a seLp£ orthogonal poly- 
nomials (Appell polynomials Ai{xi, X2, X3) for the nucleoncH13) . The remaining 
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Figure 1: Some selected models for the nucleon wave function. The model amplitudes are 
Ijolynornial approximations of second and third degree from Refs. ^(COZ), |^ (CZ), ^ (KS), 
I (GS), g (het), Q (III, IVc), and Q (1). </)as = 120x1x2x3 denotes the asymptotic form. 



A 



A 



Figure 2: A typical lowest-order contribution to virtual compton scattering in the hard- 
scattering picture of pQCD. While the hard gluon exchanges can be calculated porturbatively, 
the nucleon QDA that contains all the non-perturbative information is usually obtained 
from QCD sum-rules. 
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task was to find the corresponding expansion coefficients From the defini- 
tion of the moments of 0jv 

{ijk) := {x\x^2^^) := J [dx] x\ xi^ X3 • 4>n{xi,X2,X3) (2) 

00 „ 

= X!*^' ' / t*^^] ^1 ^2 ^3 ■ ^i(a;i,X2,a;3) , 

1=0 

which weight the QDA with different powers of the Xi, it is clear that there 
exists a linear relation between the moments of the QDA and its expansion 
coefficients Ci. ft is thus obvious that a simple matrix inversion yields the 
coefficients Ci as linear functions of the moments of 4>n- If the moments were 
known exactly for all orders, this would allow to reconstruct the QDA precisely. 

However, the non-perturbative methods presently available for actual cal- 
culations, such as QCD sum rules or QCD lattice gauge theory, can only pro- 
vide the lowest order moments up to a quite limil£d accuracy of up to 30%. 

Furthermore, as was shown in previous workald, the relation is extremely 
sensitive to uncertainties or errors in the moments. This is so because the 
matrix which has to be inverted, has a nearly vanishing determinantllll or - 
equivalently - some very small eigenvalues^. The problem is ill-posed. 

The uncertainties, which can be as large as 30%, destroy the information 
about the coefficients, e.g. by a sign flip. Therefore, the model QDAs not only 
show very different shapes for the same range of the moments input, but in 
addition the polynomial expansion does not even converge. Thus, there is an 
increase of the oscillations with increasing degree of the expansion, instead of 
an expected gradual refinement. 

So what can be done in order to circumvent the above difficulties ? 



2 A simpler approach 

In order to obtain a more reasonable model of the QDA which avoids the 
problems of unphysical oscillations (due to a failure in the fine-tuning of higher- 
order expansion coefficients) , we list some criteria of simplicity that^ physical 
distribution amplitude should fulfill in addition to QCD sum rulesLj: 

1. functional simplicity (e.g. an exponential ansatz) 

2. minimum number of parameters 

3. smooth 
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Figure 3: The haplousterotic model for the nucleon quark distribution ampHtude shown 
from different directions. This model amplitude is of exponential type — similar to a model 
used in Ref . ha - and was determined using the QCD sum-rule moments up to third order of 
Ref. h . It reaches the edge everywhere with slope 0, although this may not be evident from 
the figure. 



4. no oscillations 

5. positive 

6. substantially non-polynomial 

7. no specific process (experiment) as input (process-independence) 

A model constructed in such a way will be called "haplousterotic" (from the 
greek word aTTXovaTepo<; for "simpler"). 

Using the above criteria, our haplousterotic model amplitude 4'^'^^ was 
determined from the QCIL sum-rule moments of COZLI. The model is shown 
in Fig. |. It has the formB 

/ uir) ,{r) ,{r) \ 

^l^+ia:)=Ne^V + \ + , (3) 



with TV = 14.626, b[^^ = 0.5880, b'^^ = 8.724 x IQ-^^ and = 0.02413 for 
r = l.f| 

Of course, for the purpose of investigation of the Q^-evolution of the QDA 
or its convergence properties, ^^'^^ can be expanded into a series of Appell 
polynomials. A graphical representation of different orders of expansion is 
shown in Fig. ^ As it should be, one can observe a nice convergence. But 
one also has to notice that very high polynomial degrees are needed for the 
expansion to resemble the shape of the exact model i^^v'*^ . On the other hand, 
it is interesting to note that the second degree approximation to (f)^^~^ looks 
very much like the early model of CZ (cf. Fig. |l|). 

"There was a factor 3 missing in the b^J"^ in the first reference of Ref. 
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Figure 4: The haplousterotic model amplitude is expanded into Appell polynomials. How- 
ever, to resemble the shape of the exact (oo) form, polynomials of very high degree are 
needed. 

Finally, we would like to emphasize that a fixed model amplitude of the 
nucleon can be used to study systematically the effects of various uncertainties 
and approximations that are necessary to apply perturbative QCD in the cal- 
culation of physically relevant quantities like form factors, decay widths etc.. 
Relevant problems in this_context are e.g. constaiKl—pr dynamical pJland re- 
lated end-point problemalil, choice of Aqcd, k±-U^ and SukakovMtZl effects, 
fc^-dependence of the full wave function, higher-twist effects, higher-order as 
etc. 

3 Summary 

Due to the sensitivity of the relation between moments and coefficients to 
small variations in the moments (which are unavoidable when the moments 
are calculated, e.g. by the method of QCD sum rules) strong oscillations result 
in the polynomial approximation, which lead to a non-convergence. 

A way out of this dilemma is to restrict the structure of the quark distribu- 
tion amplitude to physical and phenomenological criteria of simplicity unless 
significantly more accurate methods for obtaining QCD wave functions are 
available. 

Our haplousterotic model of the nucleon QDA contains all available and 
physically relevant information and it will be an interesting challenge to test 
it in VCS. 
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